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In two nearby atoms, the dipole-dipole interaction can couple transitions with orthogonal dipole 
moments. This orthogonal coupling accounts for a number of interesting effects, but strongly de- 
pends on the geometry of the setup. Here, we discuss several setups of interest where the geometry 
is not fixed, such as particles in a trap or gases, by averaging over different sets of geometries. Two 
averaging methods are compared. In the first method, it is assumed that the internal electronic 
evolution is much faster than the change of geometry, whereas in the second, it is vice versa. We 
find that the orthogonal coupling typically survives even extensive averaging over different geome- 
tries, albeit with qualitatively different results for the two averaging methods. Typically, one- and 
two-dimensional averaging ranges modelling, e.g., low-dimensional gases, turn out to be the most 
promising model systems. 

PACS numbers: 42.50.Fx, 42.50.Lc, 42.50.Ct 



I. INTRODUCTION 

Two nearby atoms can interact in an energy-transfer 
process via the vacuum where one of the atoms is de- 
excited whereas the other atom is excited [H 0, S H- 
This dipole-dipole interaction has been studied in great 
detail, albeit mostly for the case of two-level atoms with 
parallel transition dipole moments 0, 0, 0, H, H, [13, [iHj 
El, [H, Q [H, d, E [3- It is known to modify the 
collective system dynamics and thus virtually all observ- 
ables considerably, as was also shown in a number of 
related experiments [H, ^ [2l|, [H, [H, [H, [H . Recently, 
it was found that a new class of effects arises from the 
dipole-dipole coupling between transitions with orthog- 
onal dipole moments [ij, [S 13 [l3]- This coupling is 
somewhat surprising since for single-atom systems, only 
near-degenerate non-orthogonal transitions can be cou- 
pled via the vacuum But in real atoms, e.g., tran- 
sitions from one state to different Zeeman-sublevels of a 
different electronic state typically have orthogonal tran- 
sition dipole moments. Therefore, the vacuum-coupling 
of such transitions in single atoms usually docs not oc- 
cur, which is unfortunate, since the corresponding cou- 
plings are known to give rise to many fascinating appli- 
cations Q. 

In contrast, orthogonal transition dipole moments in 
different atoms do interact via the vacuum, with cou- 
pling coefficients dependent on the relative alignment of 
the atoms, see Fig. [1] It was shown in IJ] that this 
interaction creates coherences involving excited states 
that are not driven by any laser fields. This observa- 
tion can be generalized by studying the two-particle mas- 
ter equation under rotations of the inter-atomic distance 
vector [l^. It was found that because of the orthogo- 



nal couplings, typically complete Zecman manifolds have 
to be considered in modelling the dipole-dipole inter- 
action of two atoms, such that the usual few-level ap- 
proximation is no longer possible. The orthogonal cou- 
plings crucially influence the system dynamics. For ex- 
ample, the long-time dynamics of a two-atom system can 
strongly depend on the relative orientation of the two 
atoms [1^. For a suitable laser and detector setup, un- 
dampcned periodic oscillations in the fluorescence inten- 
sity are observed for some relative orientations of the 
two atoms, whereas the system evolves into a station- 
ary steady state for other relative orientations. The rea- 
son for this geometry-dependence is the structure of the 
dipole-dipole constants. If the coupling of orthogonal 
transition dipole moments vanishes, then also the oscil- 
lations in the long-time limit vanish. 

In many situations of interest, however, the geometry is 
not fixed. For example, in a linear trap, the inter-atomic 
distance usually can be described classically as a sinu- 
soidal oscillation around a mean distance. In this case, 
a dependence of the dynamics on the orientation of the 
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FIG. 1: (Color online) In the left subfigure the geometry of 
our system is shown. The inter-atomic distance vector r\2 
is parameterized by the angles 6 and (f) and the length ri2. 
Atom A is located in the origin and atom B at ri2. Each 
atom is a three-level system in A configuration (b). The two 
lower states have an energy separation 5. fii {^2) is the Rabi 
frequency of the driving laser field coupling to transition 1 ^ 
3 (2 ^ 3) and the spontaneous decay rates are 71 and 72. 
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dipole moments relative to the oscillation direction can 
be expected. A gas of atoms corresponds to a setup where 
both the orientation and the distance of any given pair 
changes with time. Thus the question arises, whether the 
geometry-dependent effects of the dipole-dipole interac- 
tion of orthogonal transition dipole moments survive an 
averaging over different geometries. 



Therefore, here we discuss the fluorescence intensity 
emitted by a pair of three-level A-type atoms when av- 
eraged over sets of different geometries of interest, see 
Fig. [TJ Our primary interest is the question whether the 
dipole-dipole couplings of orthogonal transition dipole 
moments survives an averaging over different geometries 
and thus also is of relevance if the two atoms are not fixed 
in space. Since the modulations in the fluorescence inten- 
sity are a direct consequence of these couplings, they are 
a convenient indicator and allow for a quantitative anal- 
ysis. The second major question involves the way the av- 
eraging should be treated theoretically. For comparison, 
we discuss two different ansatzes. First, one can assume 
that the internal electronic dynamics is much faster than 
the change of the geometrical setup. On the other hand, 
we consider the case where the change on geometry is fast 
enough such that the atoms essentially see an averaged 
interaction potential. The latter approach for example 
is used in the context of ultracold quantum gases to de- 
rive the 1/r long-range potential from the dipole-dipole 
coupling of parallel dipole moments by averaging over 
all possible orientations of the inter-atomic distance vec- 
tors 



We find that in general the orthogonal couplings can 
survive an extensive averaging over different geometries 
as long as the inter-particle distance remains small. The 
magnitude of the effects in the averaged signal, however, 
strongly depends on the averaging range, and also on the 
averaging method. Typically, one- or two-dimensional 
systems can be expected to show larger effects of the 
dipole-dipole coupling. We also show that the two aver- 
aging methods considered can give very different results 
when averaged over the same set of geometries. In most 
situations, however, the case where the change in geome- 
try is slow as compared to the internal dynamics is more 
favorable. 



The article is organized as follows: In Sec. Ill Al we 
present the model system, derive the equations of motion 
and discuss our main observable, the time-dependent flu- 
orescence intensity. In Sec. IIIBl the two averaging meth- 
ods are presented and discussed. Sec. IIIII presents the 
results from the averaging for various different situations 
of interest. Finally, our findings are discussed and sum- 
marized in Sec. IIVI 



II. THEORY 

A. The model system 

We consider a system consisting of two identical three- 
level atoms in A configuration, see Fig. [1] The atomic 
states have energies htOi {i G {1,2,3}). The transition 
dipole moments of each individual atom are assumed 
perpendicular, as it is common for near-degenerate elec- 
tronic states in atomic systems such as Zeeman sub- 
levels. For simplicity, both transition dipole moments 
are assumed to be real; the one of the 1 <-> 3 transi- 
tion di = ((ii,0, 0)"^ is orientated along the x direc- 
tion and that of the 2 <-*■ 3 transition d2 = (0, c?2,0)^ 
along the y direction. A comparison with the case of 
complex dipole moments coupling to circularly polar- 
ized light was given in [l^ . It should be noted that it 
was found in [l6[ that in general all Zeeman sublevels 
of two nearby dipole-dipole interacting multilevel atoms 
have to be considered in order to correctly account for 
the different dipole-dipole couplings occurring in the sys- 
tem. Couplings to certain Zeeman sublevels can be elim- 
inated, however, in special geometries, or via a detuning 
between the different transition frequencies, thus recov- 
ering the well-known few-level systems. In the follow- 
ing, we are interested in arbitrary geometries, and are 
thus restricted to an elimination via detunings. A A- 
type level scheme could be realized, for example, in a 
four-level J=l/2 ^ J=l/2 scheme [2^ subject to a 
static magnetic field, such that the energy spacing be- 
tween the upper states is sufficiently large to neglect 
dipole-dipole coupling to one of the upper states in the 
four-level scheme. The frequency difference between the 
two lower states is denoted by S. Atom A is located in 
the origin of our coordinate system ri = (0,0,0)"^ and 
atom B at r2 = ri2 = ri2 (sinfl cos 0, sin0 sin cos 0)^, 
where the distance vector between the two atoms is ri2. 
The driving laser fields propagate in z direction. For this 
system the Hamiltonian reads 

(1) 
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(2b) 
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Ha represents the free energy of the atomic states. The 
free energy of the vacuum field is described by Hf. Hyac 
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is the interaction Hamiltonian of the vacuum field, and 
H]^ is the term describing the interaction with the laser 
fields in rotating- wave approximation (RWA). The laser 
fields have amplitudes £i, frequencies Vi and polariza- 
tion unit vectors (i e {1,2}), respectively, fii(r) = 
Vli exp[ifci • r] with Vli = [di ■ ei)£i/h are the correspond- 
ing Rabi frequencies. E{r) represents the quantized vac- 
uum field modes. Furthermore, ujkx is the frequency of a 
vacuum field mode with creation and annihilation oper- 
ator a^_^ and a^x , respectively. The energy of the atomic 
state \i) is fiiUi. We define atomic operators 



S 



\^)kk{j\ (i, J e {1,2, 3} and fee {1,2}), (3) 



one atom and the perpendicular one of the other atom is 
described by the expression containing the cross coupling 
constants FjJ^ and OjJ^. The symbol denotes the other 
atom than /x, e.g., for /i = 2 one has -^fi = 1. Note that 
the interaction picture in Eq. ([4]) is chosen such that the 
residual explicit time dependence exp[±iAt] which can- 
not be transformed away is attributed to the terms that 
describe dipole-dipole coupling of orthogonal transition 
dipole moments. This choice is motivated by the phys- 
ical origin of this time dependence, which arises from 
these orthogonal couplings [lB| . The frequency A is de- 
termined by A = S + A2 — Ai ~ 1^2 ~ i^i- 
The spontaneous decay rates are given by 



where \i)k denotes the ith electronic state of atom k. For 
i = j, Eq.Q corresponds to a population, whereas for 
i^j it is a transition operator. 

Choosing a suitable interaction picture it is possible 
to describe the system by the master equation for the 
atomic density operator p given by 
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and the di pole -dipole coupling constants can be calcu- 
lated from [ll] 
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In evaluating these coupling constants we have approxi- 
mated W31 (jJ32 ~ ujq. Re and Im denote the real and 
imaginary part of the tensor x whose components are 
given by 
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Sfj^i, is the Kronecker delta symbol. For our choice of the 
atomic system, the coupling constants between orthogo- 
nal dipole moments evaluate to (ry = fco ^12) 



- E (*^^c S^32^S[^^\p 



(4) 



Here, the RWA and the Born-Markov approximation 
were used. The first term, which contains the detunings 
Aj = Vi — (w3 — iOi) of the driving laser fields, appears be- 
cause of the chosen interaction picture. The interaction 
with the laser fields is expressed by the second summand 
with the Rabi frequencies Qj{r^i). The contribution con- 
taining 7j represents the individual spontaneous decay of 
each transition in the two atoms. In our case the spon- 
taneous decay rate on transition 3<->j is denoted as 2jj. 
The term with F^'' contains the dipole-dipole coupling 
between a dipole of one atom and the corresponding par- 
allel dipole of the other atom. The contribution propor- 
tional to ^I'j'^ represents the corresponding dipole-dipole 
energy shift. The interaction between a dipole moment of 
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Our main observable is the total time-dependent flu- 
orescence intensity emitted by the two atoms. It is as- 
sumed to be measured by a detector placed on the y-axis 
at the point R = RR with R = (0, 1, 0)^. This intensity 
is proportional to the normally ordered one-time correla- 
tion function 



I = {■.E'^~\R,t) E^+'>{R,t): 



(9) 
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where E^^\x,t) are the positive and negative fre- 
quency parts of the vacuum field E{x,t) = E^~\x,t) + 
E^^\x,t). For our arrangement of the detector, the 
atoms and the laser fields the fluorescence intensity re- 
duces to [la 



(10) 



where wi = (cjfj (ii)/(47reoc^-R) is a pre-factor that we 
neglect in the following. 



B. Averaging over different geometries 

The master equation Eq. ^ contains an explicit time 
dependence which is determined by the two driving laser 
field frequencies. Thus, in general, it cannot be expected 
that the system reaches a stationary steady state. This 
was demonstrated in 



15[, where it was shown that for 



A 7^ in general it depends on the relative alignment of 
the two atoms whether the system reaches a stationary 
state or not. For some geometries, the long-time limit is 
constant, whereas for other geometries a periodic oscilla- 
tion in the fluorescence intensity is predicted. Since the 
relative positions of nearby atoms in many experimental 
situations of relevance are not flxed, the question arises 
whether any time dependence survives when averaging 
over a set of geometries. The most obvious example for 
this is a three-dimensional volume of gas, where arbitrary 
relative orientations and distances can be observed. But 
also other sets of geometries may be considered. For ex- 
ample, in (23 |. an essentially one-dimensional ultracold 
quantum gas was studied. In this case, an external static 
field can be used to vary the relative alignment of dipolcs 
and the trap axis. 

In the following, we discuss two different approaches 
for calculating the averaged total fluorescence intensity, 
which is our main observable. 



1. The adiabatic case method 

In general, we have to average over the angles 0, as 
well as over the distance ri2. We discretize the respec- 
tive interval of each geometric parameter in Ni equal 
steps of size A^, respectively, where i G {r,9,(j)}. This 
gives rise to NrN^Ng different geometries. For each 
of these geometries, we evaluate the coupling constants 
and numerically integrate the master equation Eq. ([¥]). 
From this, we obtain the time-dependent fluorescence in- 
tensity llyit)]^ ^ ^ for this particular geometry (m e 

{1, . . . , A^j}). Finally, we average over all time evolutions 
of the different geometries using the expression 
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Q is a normalization constant. We work in a spherical 
coordinate system and do not only consider uniform mo- 
tions of the atoms. Thus an appropriate volume element 
^r,e,4> has to be considered. In the discretized form, the 
contributions from the different coordinates are given by 



AVr =^ Ar , 

AVe(n^) = r„^A(,, 
AV^{nr,ng) = r„^ sin(6'„JA^ 



(12a) 
(12b) 
(12c) 



When we average over one or two parameters only we 
omit the other summation(s) and volume element (s). 

This method of averaging describes the experimentally 
observable signal as long as the change of the geomet- 
ric setup is slow compared to the internal dynamics of 
the system. Then, the internal dynamics adapts to its 
long-time evolution on a timescale much faster than the 
change of the geometry. 

In the following, we will call this way of averaging the 
adiabatic case (AC) method because of the slow change 
of the geometry. 



2. The average potential method 

In our second method of averaging, a different physical 
situation is considered. Here, the change of the geome- 
try is considered fast compared to the internal dynamics. 
Then, the time evolution of the atomic system according 
to the master equation Eq. Q is not governed by cou- 
pling constants corresponding to a particular fixed geom- 
etry. Rather, the atom experiences an averaged coupling 
constant. Therefore, in this case, we start by averaging 
all coupling constants Eqs. (O over the range of geome- 
tries considered. This can be done analytically without 
a discretization of the averaging range, but again taking 
into account an appropriate volume element. Then the 
averaged coupling constants are given by 
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C dVr dVe dV^ , 



(13) 



with C e {Tf'^,nf'^,r^J,nii}. in order to average, e.g., 
over a sinusoidally oscillating distance we parameterize 
^12 = Tm + fa siu a by a mean distance r„i and an oscil- 
lation amplitude r^. In this case dVr = da, dVe = ri2d9 



and dV^ = ri2 sin 6d(j). Then, the master equation is 
solved and the fluorescence intensity is calculated using 
these averaged coupling constants. Finally, the time- 
dependent intensity is plugged into Eq. piT]) . Since the 
expression for the fluorescence intensity Eq. (fTO|) also de- 
pends on the orientation of the inter-atomic distance vec- 
tor, we also average this expression over the same set of 
geometries. 
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FIG. 2: (Color online) Time dependence of the fluorescence 
intensity averaged over ri2 by the AC method, rm = 0.25A 
and (a) Va = 0.02A, (b) = 0.I4A and (c) = 0.2A. The 
inter-atomic distance vector is oriented such that (f> = -k/A 
and 9 = 7r/2. The laser parameters are f2i =87, ^2 = 87, 
Ai = 0, A2 = 27, and the two lower states are assumed 
degenerate 5 = 0. 



FIG. 4: (Color online) Time-dependent fluorescence signal for 
different fixed distances ri2 without any averaging, (a) ri2 = 
0.10 A, (b) ri2 = 0.08 A, (c) ri2 = 0.06 A, (d) ru = 0.05 A, 
(e) ri2 = 0.04 A. The other parameters are as in Fig. [2] The 
vertical lines allow to easily judge the relative phase shifts of 
the different curves. 
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FIG. 3: (Color online) Dependence of the amplitude of the 
oscillating fluorescence intensity on the oscillation amplitude 
ri2 of the atom for (a) r,„ = 0.2A and (b) = 0.25A. The 
other parameters are as in Fig. [21 



In the following, this way of averaging will be referred 
to as averaged potential (AP) method. 



III. RESULTS 

We now turn to a numerical study of our system as 
outUned in the previous section. Different ranges of av- 
eraging will be considered, according to different setups 
of interest. In all cases, the two ways of averaging the 
fluorescence intensity will be compared. We choose as ini- 
tial condition both atoms to be in state |3) unless noted 
otherwise. 



A. Averaging over inter-particle distance 

In this section, the orientation of the inter-atomic dis- 
tance vector is fixed, while we assume a sinusoidal os- 



cillation of the distance ri2 around a mean distance 
with amplitude r^, i.e., (a) = -f sin(a) with 
a G [0,27r]. This corresponds to, e.g., atoms in a linear 
trap. In Fig. [2l we choose a mean distance = 0.25A 
and orientation 6 — 7r/2, = 7r/4. The different curves 
correspond to oscillation amplitudes 0.02A, 0.14A and 
0.2A, respectively. All curves in this figure were obtained 
using the AC method. It can be seen that in the averaged 
signal, the system does not reach a steady state in the 
long-time limit for any of these oscillation amplitudes. To 
analyze the oscillations in the long-time evolution in more 
detail, we determine the maximum (minimum) fluores- 
cence intensity /max (/min) in the long-time limit where 
the intensity undergoes periodic changes. We define an 
oscillation amplitude of the intensity as A/ = /max~/min- 
From Fig. [2l it is clear that A/ depends on the oscilla- 
tion amplitude Tq. This dependence is depicted in Fig. [3] 
for small mean distances r^, where it can be seen that 
A/ exhibits a resonance in the plot versus the oscilla- 
tion amplitude r^. This resonance can be understood as 
follows. First, one has to note that as long as the inter- 
atomic distance is not too small, typically the oscillation 
amplitude decreases with increasing particle distance, be- 
cause the coupling constants between orthogonal dipole 
moments decrease. Therefore, small inter-atomic dis- 
tances lead to a larger oscillation amplitude. Only for 
very small distances, the oscillation amplitude as well 
as the total fluorescence signal are attenuated because 
the dipole-dipole energy shifts move the atomic transi- 
tions out of resonance with the driving laser field, such 
that the upper state population is decreased. This ex- 
plains why the averaged oscillation amplitude decreases 
from the resonance maximum towards smaller oscillation 
amplitude Va- With smaller Tq, only larger inter-atomic 
distances are considered in the averaging, and thus the 
average oscillation amplitude decreases. The decrease 
of the oscillation amplitude A/ from the resonance to- 
wards higher amplitudes is due to a different mechanism. 
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FIG. 5: (Color online) Dependence of the amplitude of the 
oscillating fluorescence intensity on the oscillation amplitude 
ri2 of the atom for larger mean distance Vm = 2.25A. In 
(a) we used the AC and in (b) the AP method. The other 
parameters are as in Fig. [5] 



In Fig. [21 for both mean distances , this occurs if Va 
is large enough such that inter-atomic distances below 
about 0.06 A are included in the averaging. Some exam- 
ples of unaveraged time-dependent signals for different 
inter-atomic distances are shown in Fig. 2) For distances 
larger than about 0.06 A, the relevant contributions os- 
cillate approximately in phase, see curves (a) and (b) in 
Fig. m For smaller distances, however, the contributions 
move out of phase, as can be seen from curves (c)-(e). 
Curves (d) and (e) approximately have maxima where 
curves (a) and (b) have minima, and vice versa. Curve 
(c) is an intermediate case. Therefore, the oscillations 
with different phases cancel each other in the averaging 
process if distances below about 0.06 A are included in 
the averaging. 

In Fig.[nja) we show AI in dependence of Va for a larger 
mean distance = 2. 25 A, and over a broader range of 
oscillation amplitudes. It can be seen that the curve ex- 
hibits a series of resonances similar to the one shown 
in Fig. [21 These again occur due to an alternating de- 
structive and constructive superposition of the different 
oscillations in the averaging process. The overall ampli- 
tude AI, however, is small because of the overall larger 
inter-atomic distances considered in this figure. 

Finally, we discuss the time-averaged intensity ob- 
tained from the AP method of averaging. Some examples 
are shown in Fig. [6l Curves (a) and (c) show our results 
from the AC method and (b) and (d) those from the AP 
method. The oscillation amplitudes are ra = 0.14A and 
0.2A, respectively. All other parameters are as in Fig. [51 
The left panel shows values for 5 < j ■ t < 25 since the 
stationary oscillation is reached rapidly for these param- 
eters. Also in the AP case, the fluorescence intensity 
undergoes periodic changes in the long-time limit, see 
Fig. [6l For small oscillation amplitudes ra there is lit- 
tle difference between the two methods, see curves (a) 
and (b). However, for larger values of the amplitude 
of the oscillations in the AP case is much larger than 
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FIG. 6: (Color online) Comparison of the time dependence 
of the fluorescence averaged over r by both methods. The 
orientation angles 9 = O.Stt and <j) = 0.257r are fix and Vm = 
0.25A. (a) AC method, ra = 0.14A, (b) AP method, = 
0.14A, (c) AC method, ra = 0.2A, (d) AP method, ra = 0.2A. 
All other parameters are as in Fig. [2] 





0.22 r- 




0.2 - 








0.18 - 








0.16 - 








0.14 - 




0.12 - 




FIG. 7: Time dependence of the fluorescence averaged over 9 
by using the AC method. The distance ri2 = 0.1 A is flxed. 
(a) (j) = 0.27r, (b) (j) = Q.Stt, (c) (j) = O.Stt and (d) <p = iv. AU 
other parameters are as in Fig. [2] 



those obtained in the AC averaging. The dependence 
of the oscillation amplitude on the averaging range for 
the AP method is shown in Fig. [5l curve (b) . As in the 
corresponding curve (a) for the AC averaging method, 
resonance structures appear. But depending on ra, the 
two methods yield either similar or very different oscil- 
lation amplitudes. In addition, the result for the AP 
method seems to have a root at about ra — 0.4A. A 
careful analysis shows, however, that this minimum is 
not a true root. The reason for the minima in the AP 
curve is that for these oscillation amplitudes, the turning 
point at minimum inter-atomic distance is close to a dis- 
tance where the coupling constants between orthogonal 
transition dipole moments are small. Then, the averaged 
coupling constants are small such that the oscillation am- 
plitude has a minimum. These minima nicely show a 
crucial difference between the two averaging methods. 
In the AP method, it is easy to find averaging ranges 
where the averaged coupling constants are small or even 
vanish. Then, also the oscillation in the long-time dy- 
namics is negligible. The results from the AC method, 
however, typically remain oscillatory even for such aver- 
aging ranges, as the dynamics for each of the different 
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FIG. 8: (Color online) Time dependence of the fluorescence 
intensity averaged over 9 with fixed inter-atomic distance 
ri2 = 0.1 A. (a) AC method with cj> = 0.67r, (b) AP method 
with (j) = 0.671, (c) AC method with (j) = O.Qtt, and (d) AP 
method with <j) = O.Ott. The other parameters are as in Fig. (2] 



geometries contributes rather than only an averaged ge- 
ometry. We will find this difference again in the following 
sections. 
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FIG. 9: (Color online) Time dependence of the fluorescence 
intensity averaged over (f>. ru is fixed at 0.1 A. (a) 6 — 0.257r, 
(b) e = O.Stt, (c) e = O.Stt and (d) 6 ^ tt. Curves (a-d) 
are obtained using the AC method, (e) shows a result using 
the AP method for 6 = O.Stt. All other parameters are as in 
Fig.H 



B. Averaging over relative orientation 

In the following, we consider the case where the inter- 
atomic distance ri2 is fixed, but the relative orientation 
and thus the angles 9 and/or (j) are averaged over. A 
realization for this could be a Mexican-hat-like potential 
where one of the atoms is placed in a potential dip in the 
center whereas the other atom is confined to the poten- 
tial minimum in the rim. First we fix the angle and 
assume atom B to move around A on a circle in a plane 
which is perpendicular to the x-y-plane and includes the 
origin. Since is only defined between and tt we have 
to average over two semicircles with (j) and (f> + tt to in- 
clude the whole circle. Some examples of our results for 
ri2 = O.IA using the AC method of averaging are shown 
in Fig. [71 Here, the angle (j) is chosen as 0.27r, O.Stt, O.Stt 
and TT, respectively. For (j) = O.Stt and for = tt the 
system reaches a time-independent steady state in the 
long-time limit. That is because the cross-coupling con- 
stants are zero for these values of <f>, as both F.^J? and fJjJ^ 
are proportional to sin(20) sin^(6'), see Eqs. (|S]). But even 
though the coupHng constants are zero in both cases, the 
resulting intensities are not identical. This demonstrates 
that it is not sufficient to analyze the coupling constants 
alone to understand the system dynamics. 

In addition we can see from Fig.[7]that there is a phase 
shift of TT with respect to the oscillation in the long-time 
limit between the two curves for = 0.2tt and (f) = O.Stt. 
We found that in general curves for different values of 
(j) split into two groups separated by such a phase shift 
of TT. The first group contains curves for < (p < tt/2, 
whereas the other consists of curves for 7t/2 < <f) < tt. 
Within each of these groups, the oscillation amplitude of 
the intensity has the same dependence on the angle <f>. 
For (j) = and cj) = tt/2 the amplitude is zero. Then 
it increases with growing (f) and reaches a maximum for 
(j) = 0.2tt and 4> = O.Stt. respectively. Thus, in Fig. [7] the 



curves with maximum oscillation amplitude are shown. 
This separation is likely to appear due to the change of 
sign of the cross-coupling constants at = O.Stt. This 
can be understood from the explicit expressions of the 
couphng constants Eqs. ([5]). In the master equation 
we can see that a sign change in the terms with the cross- 
couplings can be rewritten as a constant phase shift factor 
of exp[iTr]. It is, however, not straightforward to connect 
this phase shift to the phase shift seen in Fig. [3 because 
the oscillation frequency of the time-dependent fluores- 
cence in general does not only depend on A, but also, 
e.g., on the laser field Rabi frequencies. In addition, one 
has to note that the geometric parameters 9, (j) enter the 
total fluorescence intensity as well, see Eq. (fTO)) . But 
our interpretation is further supported by the fact that a 
change of the inter-atomic distance ri2 has no influence 
on the separation of our curves into two groups. The sep- 
aration also persists for different initial conditions, e.g., 
atom A in state |1) and atom B in |3), and thus is not a 
consequence of the initial dynamics until the steady state 
has been reached. 

Using the AP method the separation into two groups 
remains, but again the curves are different from our re- 
sults from the other averaging method. In Fig. [8] we com- 
pare curves from both methods of averaging for averag- 
ing over 9 with fixed angles (j) = O.Gtt and (j) = O.Qtt. The 
curves resulting from the AP method have pronounced 
local extrema in each oscillation period in addition to the 
global ones, and the overall intensity is higher as in the 
AC case. 

Next we assume atom B to move on a circle in the 
x-y-plane. Thus 9 is fixed and we average over the an- 
gle (j). The inter-atomic distance is O.IA. Some examples 
of our results are shown in Fig. [9l where 9 is chosen as 
0.2Stt, O.Stt, O.Stt and tt. For the AC method the system 
does not reach a time-independent state in the long-time 
limit except for the angle 9 = tt. This is because for this 
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FIG. 10: (Color online) Time dependence ol the total fluo- 
rescence intensity averaged over 9 and for the inter-atomic 
distances (a) ri2 — O.IA, (b) ri2 = 0.15A and (c) ri2 = 0.2A. 
Curves (a-c) are obtained using the AC method, (d) is the 
result from the AP method for ri2 = O.IA. The other param- 
eters are as in Fig. [2] 
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FIG. 11: (Color online) Time dependence of the total fluores- 
cence intensity averaged over r, and for atom B moving on 
a sphere around atom A with additional harmonic oscillation 
of the inter-atomic distance. Here, = 0.2A and Va = 0.12A. 
In (a) we used the AP and in (b) the AC method. All other 
parameters are as in Fig. [2] 



choice of = TT the coupling constants vanish since they 
are proportional to sin(20) sin^(6'), see Eqs. For any 
different 9 our system remains oscillating in the long-time 
limit. In case oi 9 = O.Stt one can see local extrema in 
addition to the global extrema in the fluorescence inten- 
sity. In both cases, even the time-averaged average in- 
tensity is considerably larger than in the non-oscillatory 
case 9 = TT. Interestingly, for 9 = tt/2, the absolute value 
of the intensity is lower than for the non-oscillatory case 
9 = TT. Thus, the orthogonal coupling together with the 
averaging can have both an enhancing or a detrimental 
effect to the total emitted fluorescence. 

For this set of geometries, the AP method of averaging 
always yields a stationary long-time limit and thus be- 
haves qualitatively different from the first method. The 
reason for this is that the coupling constants for the or- 
thogonal couplings vanish upon averaging over then angle 
(j). As discussed before, then the time dependence in the 
long-time limit also vanishes, see Fig. [SI 

We now turn to the case of atom B moving around A on 
a sphere with radius ri2. In this case, neither of the two 
angles 9 and (j) is fixed, and we have to average over both 
of them while the inter-atomic distance is fixed. Some 
results from both methods are shown in Fig. [TOl We 
already know that the coupling constants vanish when 
averaged over 4>. That is why the time dependence in 
the AP method also vanishes when we average over 9 
and (f), see curve (d). In curves (a)-(c) obtained using the 
AC method; the inter-atomic distance is chosen as O.IA, 
0.15A and 0.2A, respectively. One can see that both the 
oscillation amplitude and the absolute value of the fluo- 
rescence intensity decrease with increasing inter-particle 
distance. For the distance 0.2A there is almost no oscil- 
lation left due to the vanishing of the coupling constants 
with increasing inter-atomic distance. This also explains 
why this curve approaches the AP method result, where 
the averaged coupling constants are zero. As compared 



to curves (a) and (b) in Fig. [SI curve (a) in Fig. [TUl shows 
that the additional averaging over 9 does lead to a re- 
duction of the oscillation amplitude. Still, the oscillation 
and thus the dipole-dipole coupling of orthogonal dipole 
moments can survive an averaging over all orientations, 
depending on the averaging case. 



C. Averaging over distance and orientation 

After the individual averaging over the inter-atomic 
distance and the relative orientation of the two atoms in 
the previous sections, we now consider the case of aver- 
aging over both. This situation is realized, e.g., in a gas 
of atoms, where the relative position of any two particles 
changes with time. An averaging over the two-particle 
configuration space is meaningful, since in a macroscopic 
volume of gas at any time there is a finite probability 
for an arbitrary geometry within the volume of the sam- 
ple to be present. A different realization is a sample of 
atoms randomly embedded in a host material. In this 
case, again an averaging is in order. The two situa- 
tions differ, however, since the former case corresponds 
to a time-dependent geometry for any two-particle sub- 
system, whereas the latter case can be represented by a 
sample of time-independent pairs. 

Thus, in the following, we investigate whether in these 
cases any time dependence of the fluorescence intensity 
remains in the long-time limit by considering a system 
where ri2, 9 and (f) are variable. The three-dimensional 
case of course leaves several possibilities for the averaging 
range. In the following, we will consider two cases. In 
the first case, atom B moves on a sphere with atom A 
in its center and additionally oscillates around the mean 
distance ri2 with an amplitude r^. In the second case, 
the particle fly-by, particle B passes atom A moving with 
constant velocity on a straight line, see Fig. [T^l 
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FIG. 12: (Color online) Geometry for the case of atom B 
flying past atom A with constant velocity on a straight line 
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FIG. 13: (Color online) Dependence of the oscillation am- 
plitude of the time-dependent intensity on Zfnax for impact 
parameter Tmin = 0.05A. (j) = -k/A, and the other parameters 
are as in Fig. [5] In (a) we used the AC and in (b) the AP 
method. 



In Sec. IIII Bl we have seen that averaging the coupling 
constants over </> makes them vanish, such that the sys- 
tem does not show any time dependence in the long-time 
limit when we use the AP method of averaging. This, 
of course also holds true for the three-dimensional aver- 
aging for atom B moving on a sphere with oscillation of 
the inter-atomic distance. In contrast, the AC method 
of averaging still yields time-dependent fluorescence in- 
tensities. An example is shown in Fig. [11] Here, the 
inter-atomic mean distance is chosen Vm = 0.2A and the 
oscillation amplitude is = 0.12A. We see that even 
if we average over all three geometric parameters, the 
system does not reach a time-independent state in the 
long-time limit, even though the oscillation amplitude is 
small. 

Finally, we consider the case where atom A flys past 
atom B along the z-axis from —Zmax to Zmax with con- 
stant velocity, see Fig [121 The angle 4> is fixed and we 
average over 9 and ri2 considering the respective vol- 
ume element. We analyzed the case (j) = n/A and found 
that for both averaging methods the fluorescence inten- 
sity remains oscillatory in the long-time limit. To further 
study these oscillations, in Fig. [121 we show the oscillation 
amplitude of the time-dependent fluorescence intensity 
in the long-time limit against the extend of the motion 
Zmax- The minimum inter-atomic distance is chosen as 
rmin = 0.05A. Curve (a) shows our results from the AP 
and (b) those from the AC method of averaging. One 
can see that in both cases the amplitude decreases with 




FIG. 14: (Color onUne) Dependence of the oscillation am- 
phtude of the time-dependent intensity on Zmax for smaller 
values of z^ax- (a) AC method with rmin = 0.05A, (b) AP 
method with rmin ~ 0.05A. All other parameters are as in 

Fig.m 



increasing z^ax for large values of Zmax- This is because 
for large distances the dipolc-dipole interaction tends to 
zero, and oscillations only occur if the particles are close. 
If the averaging interval contains increasing ranges of z 
where there essentially is no oscillation because of the 
inter-atomic distance, then the oscillations in the overall 
signal decrease. It is interesting to note, however, that 
in this averaging configuration the AP method shown 
in curve (a) yields much larger oscillations than the AC 
method shown in curve (b). Also, it can be seen that the 
AP method shows oscillations over a range of Zj^ax up to 
several wavelengths A. The reason for this is as follows. 
In the AP method, the coupling constants are averaged 
over the different geometries. For z = 0, the distance 
between the particles is rmin = 0.05A. At this position, 
the coupling constant fijj^ acquires a large value of more 
than 33O7. Of course, with increasing distance the con- 
stant fly'^ rapidly decreases down to zero. But averaging 
over a certain range [—Zmax, Zmax] still gives a consider- 
able averaged coupling constant f^.^J? even for values of 
Zmax where the unaveraged coupling constants are negli- 
gible. This is the reason why the oscillations persist for 
large Zmax values in the AP case. In contrast, in the AC 
case, contributions from larger z values do not oscillate 
at all such that the decrease of the oscillation amplitude 
with Zmax is much more rapid. 

We now focus on the region with smaller motion ex- 
tends Zmax- The corresponding results are shown in 
Fig. [H for rmi7i = 0.05A. In the hmit Zmax 0, 
the time-dependent fluorescence approaches the unaver- 
aged curves (d) and (e) in Fig. [4l which exhibit rela- 
tively low oscillation amplitudes. The reason is that at 
this small distance, the atomic states are shifted by the 
dipole-dipole interaction out of resonance with the laser 
fields, such that the overall fluorescence is low. For both 
methods, the intensity oscillations first strongly enhance 
with increasing Zmax j and then decrease again after pass- 
ing through a maximum oscillation amplitude. The AC 
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method results for larger Zmax essentially remain struc- 
tureless. The AP results, however, exhibit some oscilla- 
tions, and only then start to decay monotonously with 
increasing averaging range. Due to the complexity of the 
system, it is difficult to definitively attribute the oscilla- 
tion to a property of the system. We believe, however, 
that they arc due to a similar alternating constructive 
and destructive interference in the averaging as the one 
that led to the resonance structures in Figs.[3]and[5l Such 
resonances do not appear in the AC method results, be- 
cause there the contributions for higher values of z where 
the oscillations in the AP method appear are already too 
small. 



IV. DISCUSSION AND SUMMARY 

Dipolc-dipole interactions between transitions with or- 
thogonal transition dipole moments gives rise to a new 
class of effects in collective quantum systems. These cou- 
plings, however, strongly depend on the geometry of the 
setup, and even vanish for some geometries. Therefore 
here we have discussed different averaging schemes to 
answer the question whether measurable effects of the 
dipole-dipole coupling of orthogonal dipole moments sur- 
vive if the geometry of the system under study is not 
fixed. As observable, we chose the easily accessible fluo- 
rescence intensity of a pair of laser-driven A- type atoms, 
which for suitable laser parameters is known to exhibit 
periodic oscillations in the long-time limit due to the or- 
thogonal couplings. 

As a main result, wc found that the effects of the 



dipole-dipole coupling of orthogonal transition dipole 
moments can survive extensive averaging over all three 
spatial dimensions. We have analyzed the obtained aver- 
aged signals, and expect our physical interpretations to 
carry over to other atomic level structures. Depending 
on the averaging range, both constructive and destruc- 
tive superpositions of the contributions for the respective 
geometries is possible, such that a wide range of results 
was observed. The results also strongly depend on the 
method of averaging, and thus on the physical situation 
considered. Typically, the adiabatic case, where the ge- 
ometry changes slowly as compared to the internal dy- 
namics, is more favorable since it better preserves the in- 
tensity oscillations. In the average potential case, where 
the change of geometry is so fast that the atoms effec- 
tively see a dipole-dipole interaction averaged over the 
different geometries, some averaging ranges lead to an 
exact vanishing of the coupling constants. This usually 
does not occur in the adiabatic case. A somewhat differ- 
ent situation was found in the particle fly-by, where the 
averaging over the coupling constants in the AP method 
led to a much wider range of distances over which an 
effect of the orthogonal couplings can be observed. In 
general, our results show that the most pronounced ef- 
fects of the orthogonal couplings in systems with variable 
geometry can be expected in one- or two-dimensional se- 
tups. There, it is easier to avoid detrimental averaging 
over extended sets of geometries, and additional control 
parameters such as the orientation of the dipole moments 
with respect to the axis of a onc-dimensional sample al- 
low to study the system properties in more detail. 
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